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ABSTRACT 

We present a spherically symmetric, core-collapse model of SNR RX J17 13.7-3946 that includes a hydrody- 
namic simulation of the remnant evolution coupled to the efficient production of cosmic rays (CRs) by nonlinear 
diffusive shock acceleration (DSA). High-energy CRs that escape from the forward shock (FS) are propagated 
in surrounding dense material that simulates either a swept-up, pre-supernova shell or a nearby molecular 
cloud. The continuum emission from trapped and escaping CRs, along with the thermal X-ray emission from 
the shocked heated ISM behind the FS, integrated over the remnant, is compared against broadband observa- 
tions. Our results show conclusively that, overall, the GeV-TeV emission is dominated by inverse-Compton 
from CR electrons if the supernova is isolated regardless of its type, i.e., not interacting with a ;§> 100 M shell 
or cloud. If the SNR is interacting with a much larger mass > 10 M Q , pion-decay from the escaping CRs may 
dominate the TeV emission, although a precise fit at high energy will depend on the still uncertain details of 
how the highest energy CRs are accelerated by, and escape from, the FS. Based on morphological and other 
constraints, we consider the 10 4 M Q pion-decay scenario highly unlikely for SNR RX J1713. 7-3946 regard- 
less of the details of CR escape. Importantly, even though CR electrons dominate the GeV-TeV emission, the 
efficient production of CR ions is an essential part of our leptonic model. 

Subject headings: acceleration of particles, shock waves, ISM: cosmic rays, ISM: supernova remnants, mag- 
netic fields, turbulence 



1. INTRODUCTION 

The direct observation of radio, non-thermal X-ray, and 
GeV-TeV emission from several young supernova remnants 
(SNRs) provides unambiguous proof that these objects pro- 
duce relativistic particles. The morphology of this emission, 
which is often seen in thin, rim-like structures associated with 
the supernova (SN) blast wave, provides convincing evidence 
that diffusive shock acceleration (DSA) is the mechanism pro- 
ducing these particles. This and other evidence, such as the 
direct measure of efficiency at the Earth bow shock (e.g., 
lEllison et al.ll990h . and the likely existence of amplified mag- 
netic fields in some SNRs, suggests that the efficiency for pro- 
ducing superthermal ions, i.e., cosmic rays (CRs), is high in 
DSA. A critical issue for DSA, and for the origin of CRs, 
however, concerns the radiation mechanism responsible for 
the GeV-TeV emission. Is it dominated by pion-decay emis- 
sion from ions or inverse-Compton emission from leptons? 

Given relativistic electrons and ions, three mechanisms 
- non-thermal bremsstrahlung, inverse-Compton (IC) scat- 
tering, and pion-decay - can produce GeV-TeV photons. 
Sorting out this mix in a particular remnant not only pro- 
vides information on how CRs are produced, it constrains 
the basic physics of DSA, in particular the electron to 
proton ratio, K cp , of accelerated particles. While most 
modelers of young remnants claim that shock accelerated 
proto ns are mainly responsible for the GeV-TeV e mission 
(e.g., iMorlino et all 120091 iBerezhko & Volld l2010h some 



researchers (e.g.. iPorter et al.l l2006t iKatz & Waxmanl 120081: 
IZirakashvili & Aharonianl2010h . have suggested that IC from 
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relativistic electrons dominates the GeV-TeV emission from 
SNR RX J1713.7-3946 (henceforth SNR J1713). Recently 
we have definitively shown, with a model that includes 
the production of thermal X-ray lines s elf-consistently wit h 
the broadband continuum emission (i.e., lEllison et alJEOlOh . 
that IC emission from relativistic electrons clearly dominates 
pion-decay if SNR J 17 13 is in a homogeneous environment 
typical of a Type la SN. 

The assumption of a homogeneous environment, although 
typically assumed in self-consistent models of SNR 11713, 
is an important caveat because SNR J1713 is more likely 
a core-collapse S N that may have e xploded in a complex 
environme nt (e.g.. iFukui et alJ l2003h iMoriguchi et alJ 120051 : 
llnoue et al . 2009; S ano et al.ll201QF Because of this, we have 
generalized our model in two important ways. First, we can 
now model a spherically symmetric, non-homogeneous envi- 
ronment. We consider a SN that explodes in a pre-SN wind 
surrounded by a dense, swept-up shell of wind material. We 
further allow for external dense material that might be as- 
sociated with molecular clouds out of which massive stars 
form. Second, in addition to 7-rays produced by CRs trapped 
within the forward shock (FS), we now include 7-ray produc- 
tion from the CR protons that escape from the FS and diffuse 
rapidly into the dense shell and/or the molecular cloud mate- 
rial beyond the forward shock. 

These generalizations allow a fairly realistic model of 
a core-collapse SNR J1713 interacting with dense mate- 
rial, at the expense of adding model parameters. How- 
ever, despite the added parameters and the flexibilty in fit- 
ting they allow, we find it is impossible to produce a good 
fit to the broadband data, including the X-ray emission 
lines, with pion-decay dominating inverse-Compton at GeV- 
TeV energies unless the SNR is interacting with a dense 
mass concentration (> 10 4 M^) such as a molecular cloud 
(e.g. JGabici & Aharonianll2007t iGabici et al.ll2009h . An iso- 
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lated SNR, regardless of its type, even allowing for es- 
caping CRs interacting with a dense shell of material with 
M s he\\ ^ 100 M Q from a pre-SN wind , is exc luded. 

The conclusion of lEllison et al.l d2010l) . that inverse- 
Compton from electrons dominated the GeV-TeV emission, 
hinged on the assumption that both non-thermal pion-decay 
and thermal X-ray line emission scale in approximately the 
same way with the ambient density. This is expected to be 
the case if the shock accelerated particles are drawn from the 
target material, i.e., if the material swept-up by the FS pro- 
vides seed particles for DSA and simultaneously provides a 
dense target for proton-proton collisions. If the environment 
is homogeneous and has a uni form density everyw here out- 
side of the FS, as assumed in lEllison et al.l ( 120101) . photons 
produced outside of the FS, where the density is considerably 
less than behind the shock, can be ignored. In this case, the 
parameters chosen for the injection of electrons and protons 
into DSA, along with the ambient density and magnetic field, 
largely determine the IC/pion-decay ratio and t he thermal X- 
ray em ission/pion-decay ratio. As shown by lEllison et al.l 
(2010), reasonable parameters for SNR J 17 1 3 exclude pion- 
decay as the dominant radiation mechanism for the GeV-TeV 
emission in this scenario. 

There are three main populations of shock accelerated CRs 
that are important for producing 7-rays: relativistic electrons 
producing 7-rays through inverse-Compton and non-thermal 
bremsstrahlung; CR ions that remain trapped within the for- 
ward shock; and CR ions that are accelerated by the for- 
ward shock but escape upstreamQ These three populations 
are produced simultaneously by DSA in the blast wave and, 
while they have different properties, their radiation signa- 
tures, particularly IC and pion-decay, are not easily distin- 
guished at GeV-TeV energies with current observations. The 
recent Fermi-LAT observations of SNR J 17 1 3, however, now 
clearly support a leptonic model for the GeV-TeV emission 
dAbdo et al.ll20Tlb . 

An important distinction between trapped and escaping CR 
ions is that, in a non-homogeneous environment, where re- 
gions outside of the FS can have an elevated density, the 
escaping CR ions can diffuse into these dense regions and 
produce 7-rays without producing any corresponding thermal 
X-ray or IC emissionQ The production of these external 7- 
rays depends mainly on the amount of external mass but the 
diffusion properties of the CRs in the circumstell a r med ium 
(CSM) also come into play. In lEllison & Bykovl (120111) we 
presented a simple Monte Carlo model for escaping CR dif- 
fusion and investigated how the 7-ray flux depends on the CR 
diffusion without application to a particular SNR. For rea- 
sonable values of the diffusion, it is clear that, if the escap- 
ing CRs have enough external mass to interact with, the 7- 
ray flux from pion-decay can overwhelm t hat from IC (e.g., 
iGabici & Aharoniaril2"007t iLee et al.ll2008h . Other observa- 
tional considerations, such as the existence of nearby molec- 
ular clouds and the spatial coincidence of radio, X-ray syn- 
chrotron, and 7-ray emission must be used to discriminate 

4 A fourth particle population that we don't consider here are secondary 
electro n-positron pairs pr oduced by proton-proton interactions (see, for ex- 
ample, Gabici et al. 2009). These leptons will produce IC emission and may 
be important depending on the external mass concentration. 

5 For the purposes of this paper, we assume that only CR ions escape 
upstream from the shock. Relativistic electrons will suffer radiation losses 
and are not as likely to be an important contributor to the IC emission beyond 
the FS. Including escaping CR electrons would increase the IC/pion-decay 
ratio from the values we calculate and strengthen our basic conclusions. 



between a hadronic or leptonic origin in this case. 

Here, we model the broadband emission from SNR J 17 13 
for four non-homogenous environments. The first is for an 
isolated core-collapse SNR expanding in a slow pre-SN wind 
with an external, dense shell of swept-up wind material. For 
an isolated SNR, the mass of the swept-up wind must be 
less than the mass of the main sequence star and we use 
Afsheii = 100 M Q as an extreme case. In the second example, 
we increase the external mass to M s h. e ii = 10 4 M Q to minmic 
a SNR which is near a mass concentration but has not yet im- 
pacted it. In the third example, we allow the SNR to impact 
a dense shell of 10 4 M Q and follow the evolution of the rem- 
nant and the photon production as the FS moves into the shell. 
Finally, we consider the case of a fast pre-SN wind. 

We find, largely independent of diffusion parameters, that 
escaping CRs interacting with an external mass shell from a 
pre-SN wind with mass < 100 M Q will not produce enough 
pion-decay emission to compete with the IC emission from 
trapped CR electrons. If larger external masses > 10 4 M are 
allowed, as would be the case for a SNR near or interacting 
with a molecular cloud, the pion-decay can dominate the IC 
at TeV energies. Now, however, other considerations come 
into play such as the shape of the escaping CR distribution 
and the details of how CRs can stream into a molecular cloud 
and produce self-generated turbulence. We expect the escap- 
ing CR distribution to be far too narrow to allow a good fit 
to the GeV-TeV observations without a substantial contribu- 
tion from IC photons. Nevertheless, the shape is uncertain 
since escaping CRs will be highly anisotropic and the plasma 
physics of these distributions is an active area of research^ 

The broadband fit we obtain for our IC core-collapse 
model is, in fact, significantly better than that obtained in 
lEllison et al.l (120101) . particularly for the highest energy HESS 
points. The reason, as we show below, stems from the fact that 
the pre-SN wind magnetic field, in which the SN explodes, 
can be considerably lower than 3 /iG, allowing electrons to be 
accelerated to higher energies before radiation losses domi- 
nate. 

It is important to emphasize that our model, and all self- 
consistent shock acceleration models of SNR J 17 13 we are 
aware of, accelerate CR ions efficiently. The models we 
show below place 25 to 50% of the FS ram kinetic energy 
flux into relativistic ions at any instant. Only 0.25% or less of 
the instantaneous ram kinetic energy flux goes into relativis- 
tic electrons. Leptons dominate the observed emission sim- 
ply because leptons radiate far more efficiently than ions, not 
because ions are missing. Furthermore, our best-fit param- 
eters for this remnant result in maximum proton energies of 
~ 10 14 eV. Iron nuclei would be accelerated to ~ 26xl0 14 eV, 
well into the CR "knee." We make no claim that the spec- 
tral shape produced in this particular SNR matches the CR 
spect rum observed at Earth. As has been known for decades 
(e.g.. iGinzburg & Svrovat skii 1964), CRs observed at Earth 

6 We note that in addition to the direct escape of CRs upstream from the 
shock that we consider, an additional factor comes into play in an expand- 
ing SNR. As the remnant expands, the precursor region beyond the forward 
shock that is filled with CRs expands producing a "dilution" of CR energy 
densi ty. This effect has bee n studied i n detail by B erezhko and co-workers 
(e.g., Berezhko et al. 1996a b) (see also Drurv 2010). In a real shock, the di- 
lution effect is coupled to escape since the lowering of the CR energy density 
results in less efficient generation of magnetic turbulence and this will change 
the escape of CRs. Both the flow of energy out of the shock by escape and 
the dilution of the CR energy density influence the Rankine-Hugoniot con- 
servation relations in similar ways. Both act as energy sinks and both result 
in an increase in the shock compression and other nonlinear effects. 
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come from a mix of many SNRs, w ith some admixture of 
other sources (see IMever etall 11997b . after a long and tor- 
tuous propagation in the ISM. The spectrum produced by 
any one source will be modifed substantial ly by propaga- 
tion (see figure 1 in En son &Eichlerf[l985l) . and different 
elements, e.g., hydrogen and he lium, may propagate differ- 
ently (e.g.jBlasi & Amatoll201 lb . Despite some recent claims 
dAdriani et al.1201 11) , webelieve our detailed fit to SNR J1713 
is fully consistent with SNRs being the primary source of 
galactic CR ions at least to energies into the CR knee. 

2. MODEL 

The cosmic -ray, hydrodynamic, non-equilibrium ionization 
(CR-hydro-NEI) model of an evolving SNR we use has been 
described in detail in a number of pr evious papers . It is e ssen- 
tially the same as that described in Ellison et al.l ( 120101) and 
references therein except that we have added the escape of 
high-energy CRs from the shock precursor and the propaga- 
tion and 7-ray production of these CRs as they diffuse in the 
circumstellar medium. Our method de scribing the escaping 
CRs is given in lEllison & Bvkovl (1201 ll) . 

The important features of the CR-hydro-NEI model are: 
(i) the hydrodynamics and evolution of the SNR are coupled 
to the efficient production of CRs via nonlinear DSA; (ii) 
the NL CR acceler ation is calculated using the semi-analytic 
technique given in iBlasi et aTl (120051) : (iii) the CRs that re- 
mained trapped within the FS have a normalization and max- 
imum energy cutoff consistent with the normalization and 
shape of the high-energy CRs that escape upstream; (iv) mag- 
netic field amplification (MFA), presummably produced by 
nonlinear DSA at the FS, is simp l y para meterized; (v) as 
described in detail in lEllison et al.l (120071) . the NL acceler- 
ation model determines the heating of the shocked plasma; 
(vi) the electron temperature and non-equilibrium ionization 
state of the shocked plasma are calculated and followed self- 
consistently with the SNR dynamics in the interaction region 
between the contact discontinuity and the FS as the SNR 
evolves; (vii) within the SNR, thermal X-ray emission lines 
and continuum emission from synchrotron, inverse-Compton, 
bremsstrahlung, and pion-decay from trapped CR electrons 
and protons are calculated taking into account the SNR evo- 
lution and adiabatic and radiation lossesQ and (viii) the pion- 
decay emission from escaping CR protons and trapped CRs 
is absolutely normalized relative to each other and relative to 
the continuum and line emission in all other bands. 

Since, as in our previous applications of the CR-Hydro- 
NEI model, we assume spherical symmetry, all of our escap- 
ing CRs interact with the external shell. In 3-D geometry, 
where the external mass is concentrated in clumps, only some 
fraction of the escaping CRs would e ncounter the clumps be- 
fore d iffusing away (see, for example. iLee. Kamae. & Ellison! 
l2008h . Our 7-ray fluxes from escaping CRs interacting with a 
given external mass are therefore upper limits. 

Any realistic description of a SNR requires a large number 
of parameters. Parameters are required to model the super- 
nova explosion and SNR, the DSA mechanism, the CSM, and 
the diffusion of escaping CRs in the CSM. In Tables [T]|2j and 
[3] we list these input parameters, and in Table |4] we list some 
output values. 

7 As in previous applications of our CR-hydro-NEI model, we only cal- 
culate X-ray emission from the shocked ISM. If thermal emission from the 
shocked ejecta material was included, the X-ray line emission we show would 
be greater and our basic conclusion concerning the dominance of IC emission 
over pion-decay at GeV-TeV energies would be strengthened. 



Besides the assumed age and distance to SNR J 17 1 3, Ta- 
ble Q] gives values for the SN explosion energy i?sN> and 
ejecta mass M e j, as well as the power-law index of the initial 
ejecta mass distribution, n (see lEllison et al.ir2004 for a full 
discussion of n). All of our examples here assume a pre-SN 
wind with a speed V^nd, mass-loss rate dM jdt, and tempera- 
ture Twind, all of which are assumed constant. The parameter 
Cwind shown in Table Q] determines the wind magnetic field 
at a radius R from the explosion according to 



B (R) = 



wind 

V^dM/dt) 1 ' 2 



R 



(1) 



as in | Chevalier & Luol d 19941) and lEllison & Cassam-Chenaj 
(120051) . The constant parameter cr w i n d is the ratio of magnetic 
field energy density to kinetic energy density in the wind and 
can b e related to properties of the star by (e.g., Walder et al. 
l20ll 

BlRl fv r \ 2 
{dM/dt)v w \v w J 

Here, -B* denotes the surface magnetic field of the star, i?* 
is the stellar radius, v r is the star rotation velocity, and 
v w = V wul d is the terminal speed of the wind. Obtain- 
ing CTwind ~ 0.03 by fitting the spectrum of SNR J 17 1 3 as 
we do, constrains the p rogenitor star parameters. As noted 
by I Walder et al.l (1201 ll) . values of <7 W j n( j <C 1 indicate that 
the stellar wind dominates the magnetic field, producing a 
roughly radial field far from the star. 

All of t he parameters in Ta ble |2l have been described previ- 
ously (see lEllison et al.l2010l and references therein). Briefly, 
£dsa is the fraction of ram kinetic energy flux that is placed 
in superthermal particles at any instant^ -B a mp is an ampli- 
fication factor for the shocked magnetic field, K cp is the ra- 
tio of electrons to protons at relativistic energies, / s k is the 
fraction of the FS radius used to determine the maximum CR 
momentum, a cut is a factor that determines the shape of the 
particle cutoff at high momentum for trapped CRs, and iV cu t 
(along with a cu t) determines the shap e of the escaping CR 
distribution (see Ellison & Bvk ov 201 1, for a detailed discus- 
sion of a cu t and N cut ). The parameter / cq = 1 indicates 
that Coulomb equilibration is used for electron heating while 
/ cq = indicates that electrons equil ibrate instantly with pro- 
tons downstream from the shock (see lEllison et all2007l) . and 
/norm is the overall normalization applied to the photon emis- 
sion to match the observations. In all cases, a solar elemental 
composition has been assumed for the CSM. 

The values given in Table[3]determine the characteristics of 
th e CSM and the diffusion of the escaping CRs. As described 
in lEllison & Bvkovl (1201 ll) . the diffusive mean free path for 
the escaping CRs in the CSM is parameterized by 



Acsm = AcsM,o(' , g /?' 9 ,o) Qrs ["csM(-R)/^o 



(3) 



Here, r g = pc/(eB) is the gyroradius, ncsM(-R) is the CSM 
proton number density^ and a rg and (3 n are parameters. For 
scaling, we use no = 1cm -3 , r 9:0 = 10GcV/(ei?cSM,o)> 
and £>csm,o = 3 /iG. We note that instead of the jicsm term, 



8 W e note that we determine the acceleration efficiency in the [Blasi et all 
120051) semi-analytic calculation by setting £dsa ar, d then calculating the 
injection fraction Xini rather than the reverse, as is typically done (see 
Ellison et al. 2010). With a constant £dsa> Xinj varies over the lifetime of 
the SNR. 

9 The density ncsM is the value of n p shown in Figs. [T] f5] |5] and[7jat a 
given R outside of the FS. 
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FIG. 1 . — The top two panels show the proton number density and the mag- 
netic field as a function of distance from the center of the SNR for Model A. 
In each of these two panels, the dashed curve is the profile at the beginning 
of the simulation and the solid curve is the profile at tgNR = 1630 yr. The 
third panel shows the mass within R at t = and the fourth panel shows the 
escaping CR number density. The diffusion parameters as defined in Eq. [5] 
are listed in the fourth panel. Escaping CRs are only followed beyond the FS 
and they leave the spherically symmetric simulation freely at the outer radius 
of ~ 16 pc. The sharp dropoff in Acsm within ~ 9pc indicated in the bot- 
tom panel shows the effect of assuming Bohm diffusion for the trapped CRs. 



we could have scaled Acsm with (B / Bcsm.o)"' 3 " , or even 
a combination of density and magnetic field terms. These are 
essentially equivalent parameterizations unless the connection 
between background field, ambient density, and wave gener- 
ation by streaming CRs is specified. The normalization of 
the CSM diffusion coefficient, -Dcsm.o = Acsm.oc/3, can 
be estimated from CR propagation stud ies (see, for example, 
iPtuskin et al.l 12006; Gabici et al. 2009). For example, with 
-Dcsm.o = 10 27 cm 2 s~\ ncsM = 0.01 cm -3 , a rg = 0.5, 
an d n = 1, Aqsm ~ 1 pc at 1 GeV, consistent with the fits 
of IPtuskin et al.l (120061) . In general, the stronger the diffusion 
(i.e., the smaller Acsm) the greater the 7-ray emission will be 
in the external material. 

The values n un j and n s h e ii are the proton number densities 
for the uniform CSM beyond the dense shell and for the dense 
shell, respectively. The values M s h n and i? s hcii are the mass 
of the dense shell and its inner radius, respectively, and -B s hcii 
is the magnetic field in the shell. As shown in Figs. Q] and [3] 
we smooth the transition from the pre-SN wind to the dense 
shell. 

We note that within the FS we assume Bohm diffusion for 
the CRs with a mean free path A ~ r g which is very much 




FIG. 2.— Model A fit to SNR J1713 observations. The different emis- 
sion processes are: synchrotron (solid blue curve), IC (dot-dashed pur- 
ple curve), pion-decay from trapped CRs (dashed red curve), pion-decay 
from escaping CRs ( dotted black curve), and thermal X-rays (solid black 
curve). The dashe d black curve is the summe d emission. The data is from 
IAcero etalJ <2009) (radio ).ITanaka et a l. (2008) (Suzakii X-ravs), rAbdo et all 
1201 II) (Fermi-LAT), and lAharonian et all )201ll) (HESS). Note that the two 
lowest energy Fermi-LAT points are upper limits. For all models we use a 
column density of nu = 7.9 X 10 21 cm -2 . 

smaller than Acsm- This is reflected in the sharp drop in A 
within the FS as shown in the bottom panels of Figs.Q]and[3] 

3. RESULTS 
3.1. Pre-SN Wind Interaction 

For our core-collapse model A, we take the SN explosion 
energy to be £?sn = 10 51 erg, the ejecta mass M e j = 3 M Q , 
and assume a slow, dense, pre-SN wind with a mass-loss rate 
dM/dt — lO~ 5 M yr~\ and wind speed V win d = 20km 
s _1 . Our model is a simplified description that might resem- 
ble what happens after an early-type star with a fast wind 
creates a large, low-density bubble before evolving into a 
red-supergiant w ith a much slower wind (see, for example, 
IChevalierlll99"9h . As we show below, the critical conditions 
that result in a good leptonic fit are that the density is rela- 
tively low and the B-field in the wind is lower than the normal 
ISM field. Other than this, none of our conclusions depend 
critically on particular wind parameters. 

To determine the unshocked magnetic field as a function 
of radius, R, in the pre-SN wind, we take cr W j n( j = 0.03 
in equation (fl}. At the assumed age of SNR J 17 1 3 (i.e., 
tSNR — 1630 yr), the FS has not yet reached the dense ma- 
terial of the swept-up wind. The situation is shown in Fig.[T] 
where, in the top two panels, the proton number density, n p , 
and the magnetic field, B, are plotted as functions of radius, 
R, from the center of the SNR. The dashed curves in the top 
two panels are the density and magnetic field profiles at the 
start of the simulation. The solid curves are these profiles at 
tsNR = 1630 yr. Parameters have been chosen so the SNR 
radius is ~ 9pc at <snr = 1630 yr, consistent with a distance 
to SNR J 17 13 of ~ 1 kpc and an angular size of ~ 60 arcmin. 
As seen in the second panel, the pre-SN wind magnetic field 
just upstream of the FS, as determined by cr W i n< j, is ~ 0.2 /iG 
at isNR = 1630 yr and this is increased by compression and 
amplification to ~ 10 /iG immediately downstream. 

At a radius beyond the FS, we have placed a dense shell 
with a total mass ~ 100 M Q and the third panel in Fig.Q]gives 
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the mass within R. A shell mass of 100 M© is an extreme case 
and any shell from swept-up, pre-SN wind material would be 
considerably less. The fourth panel in Fig.Q]shows the num- 
ber density of escaping CRs that have diffused beyond the FS 
and the bottom panel shows the scattering mean free path be- 
yond the FS, Acsm, for lOGeV (solid black curve) and 1 TeV 
(dashed red curve) CRs. Note that we have set B = 3 /j,G 
in the shell and beyond. Since we don't calculate the syn- 
chrotron emission from electrons beyond the FS, the value of 
B beyond the FS is unimportant. 

In Fig. |2] we show the radiation produced by the trapped 
and escap ing CRs, along with the observations of SNR J 1 7 1 3 
at rad io (lAcero et al.l 120091). X-ray (Suz aku: iTanaka et all 
120081) . GeV (Fermi-LAT lAbdo et al.l201 ll) . and TeV energies 
(HESS ; lAharonian et al.ll2007l) . As described in lEllison et all 
i20 101) . this is a "best fit" result where we have varied pa- 
rameters to obtain a good match to the broadband data. The 
essential features of the comparison are: (i) inverse-Compton 
emission from relativistic electrons interacting with the cos- 
mic microwave background (CMB) strongly dominates the 
GeV- TeV emission; (ii) the quality of this broadband fit with a 
pre-SN wind, including weak X-ray emission lines consistent 
with the Suzaku observations, is superior t o that obtained 
with a uniform CSM model by lEllison et al.l (120101) ; (iii) the 
escaping CRs (black dotted curve), even though they diffuse 
through a region with 100 M Q of material, give a small contri- 
bution to the TeV pion-decay emission, and; (iv) even though 
there are a lot of relatively free parameters for this fit, the 
critical result that IC dominates over pion-decay at GeV- TeV 
energies is robust. 

We don't show a fit where parameters have been optimized 
for pion-decay dominance at GeV- TeV energies because the 
result is essentially the same as shown in lEllison et al.1 (1201 Oh 
for the homogeneous case. As is clear from Fig. [2] if K cp 
is decreased to allow pion-decay to dominate IC, the over- 
all density will have to be increased substantially to produce 
enough flux at GeV- TeV energies. Any increase in density 
(here, for the pre-SN wind through varying dM/dt and/or 
Vwind) will increase the X-ray line emission along with the 
pion-decay emission, producing a conflict with the Suzaku 
observations. In fact, as long as the contribution from the es- 
caping CRs remains well below the IC, all fitting parameters, 
such as the efficiency for PSA, <SV)SA, are co nstrained here 
just as in the fits shown in lEllison et al.l(l2010h . 

We note that the parameter / s k determines the precursor 
length in our model but that we do not explicitly model the 
spatial properties of the precur sor in our implemen tation of 
the Blasi et al. DSA model (see lCaprioli et alj|2009l and ref- 
erences therein, for an updated version of the model that does 
include the precursor explicitly). We do, however, account 
for all CRs that are accelerated and, except for the escaping 
CRs, assume they are all trapped behind the FS. We, there- 
fore, overestimate the pion-decay flux since the trapped CRs 
interact with the dense shocked plasma rather than the thin- 
ner precursor material. The only possible case where the pre- 
cursor CRs might enhance the pion-decay emission is if the 
precursor is interacting with external material denser than the 
shocked plasma. We consider this case unlikely because it re- 
quires fine tuning. If the external material is fully outside of 
the prec ursor, it is contained in the case we discuss in Sec- 
tion [3T2] below. If the FS is also interacting wit h the exter- 
nal material, it is the case we discuss in Section 13.31 below. 
Only if the precursor, but not the FS, is impacting dense ma- 
terial is there a possibility that the pion-decay emission will be 
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FIG. 3. — Same format as in Fig. [T]for Model B, where the mass of the 
external shell is ~ 10 4 Mq. 

greater than we estimate. With a precursor length determined 
by / s k = 0.1, we consider this to be unlikely. While larger 
values of / s k are possible, they imply strong self-generated 
turbulence far upstream where the density of accelerated CRs 
has dropped significantly from that at the shock (this is, in ef- 
fect, CR "dilution" as discussed by iBerezhko et al.lll996alfbT) . 
It is noteworthy that neglecting the spatial properties of the 
precursor is less likely to overestimate the electron contri- 
bution to the GeV- TeV emission since radiation losses will 
prevent the highest energy electrons from streaming far up- 
stream. 

There are two important improvement s over t he broadband 
IC fit given in Figure 4 of lEllison et al.l (120101) . One is that 
the newer Fermi-LAT data for SNR J 17 13 now clearly favor 
a IC m odel, whereas the pr eliminary Fermi-LAT data avail- 
able for lEllison etail (120101) were less clear. The second im- 
provement is in our match to the hig hest energy HESS po ints. 
In the constant ISM model used in lEllison et al.l d2010l) . the 
IC fit fell below the highest energy HESS points. Now, with 
our core-collapse model, we are able to fit the highest energy 
points successfully with only the CMB photon field. This 
improvement comes about because the magnetic field at the 
FS is lower in the core-collapse case and electrons can obtain 
a higher energy before synchrotron losses dominate. Mag- 
netic field amplification is still important (we fit the data with 
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amp 



8.5 for this case) but starting with a lower ambient 



field is advantageous. 



3.2. External Molecular Cloud Interaction 
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FIG. 4. — Same format as in Fig.|2]for Model B. As in Fig.H the black 
dashed curve is the total emission and, in this case, it lies above the HESS 
data. While a better fit could be obtained by reducing the external mass and/or 
by increasing the CR diffusion coefficient in the CSM to reduce the contri- 
bution from escaping CRs, a good fit to the GeV-TeV emission with only 
pion-decay isn't possible. 



The contribution from the escaping CRs to the TeV 
emission will increase as the external target material in- 
creases, as would be the case if the escaping CRs from 
the SNR diffused i n to a n earby molecular cloud (e.g. , 
Aharonian & Atovanl 119961: iPtuskin & Zirakashvilil 120051: 
Gabici & Aharonianl 120071: iGabici et al.ll2009t iCaprioli et al.l 
20 1 Ot lOhira et alj|20 1 lb . In Figs.|3]and|4] we show an exam- 
ple (Model B) where the mass external to the FS is 10 4 M© 
with a density of ~ 100 cm~ 3 . All other parameters are the 
same as for Model A and again we note that our model is 
spherically symmetric so all escaping CRs interact with the 
outer shell. Now, the pion-decay emission from the escaping 
CRs at TeV energies is well above the pion-decay from the 
trapped CR protons and is comparable at the highest energies 
to the IC from the trapped CR electrons. 

While a larger external mass would clearly result in pion- 
decay dominating the IC, it remains to be seen if a satisfac- 
tory fit for SNR J 17 13 can be found with just pion-decay. The 
first problem concerns the shape of the pion-decay emission 
from the escaping CR distribution. The distributions shown 
in Figs. |2] and [4] are much too narrow to produce a good 
fit to both the Fermi-LAT and HESS fluxes and a substan- 
tial contribution from IC is required to produce a goo d fit. 
We note, for example , that lZira kashvili & Ptuskinl d2008l) and 
ICaprioli et al.l d2010l) find similar narrow escaping CR distri- 
butions, as does IVladimirov et al.l ( 12006b with Monte Carlo 
simulations that dir ectly determine the escap ing flux, but that 
IGabici et all (120091) and lOhira et al.l (12011b assume broader 
distributions when averaged over the age o f the remnantFI 
As was discussed in lEllison & Bvkovl(l201 lb . the shape of the 
escaping CR particle distribution depends on the details of 
how the highest energy particles, trapped and escaping, gener- 
ate turbulence, particularly via long-wavelength effects. The 
long-wavelength, wave-particle interactions that are impor- 
tant for these high-energy, escaping particles involve highly 
anisotropic distributions that are only beginning to be mod- 

10 There are data from spacecraft observations of the Earth bow shock 
supp orting the direct escape of a narrow distri bution of accelerated particles 
fe.g., IScholer"e7aIi 1980: Mitchell et aU1983l) . 
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FIG. 5.— Impact of the FS on a 10 4 M spherical shell for Model C. The 
left panels show the density profiles and the right panels show the high-energy 
emission at the ages indicated. As in Figs. [T] and [5] the dashed black curves 
in the left panels show the density profiles at the start of the simulation and 
the solid red curves show the profiles at the ages indicated. For the right-hand 
panels, the various emission processes have the same line characteristics as 
in Figs. f2] and [4] i.e., synchrotron (solid blue curve), IC (dot-dashed purple 
curve), pion-decay from trapped CRs (dashed red curve), pion-decay from 
escaping CRs ( dotted black curve), and thermal X-rays (solid black curve). 
In the 1600yr right-hand panel we show the high-energy data along with 
the summed emission (dashed black curve). While the parameters for the 
1600yr case give a reasonable fit to the broadband, high-energy data with 
pion-decay (from combined trapped and escaping CRs) dominating the GeV- 
TeV emission, a better fit can only be obtained by increasing the IC relative 
to the pion-decay. The thermal X-rays are calculated assuming Coulomb 
heating of electrons. 

eled self-consistently (e.g.. ICaprioli. Amato. & Blasil I2010t 
iBykov. Osipov. & EllisorJl201 U ISchure & Belli I20H1) . It is 
certainly possible that broader distributions of escaping CRs 
than we show may come from these calculations or others, and 
that a better fit to the GeV-TeV flux without the IC component 
could be obtained. Nevertheless, we feel it is unlikely that es- 
caping CRs streaming away from a relatively young SNR will 
have a distribution much broader then we show here. 

Regardless of the details of the wave generation, however, 
the shapes of the trapped CR ions and the escaping ones 
are related since the turnover in the trapped ion spectrum 
comes about as CRs escape from the FS. In the case of a 
high enough £>-field, the shape of the electron turnover is 
determined mainly by radiation losses. With any attempt to 
broaden the spectral cutoff to allow a match to the GeV-TeV 
emission, care must be taken so that the X-ray synchrotron 
isn't broadened beyond acceptable limits. This is particularly 
important for a core-collapse model where the effective mag- 
netic field is less than in a homogeneous model. 

Another problem that must be addressed if 7-ray produc- 
tion in external material is important concerns the observed 
spatial coincidence of 7-ray and X-ray emission. If 7-rays 
from pion production in external clouds were to dominate 
the 7-ray emission one would not expect the observed good 
agreement between th e X-ray and 7-ray morphologies (i.e., 
lAharonian et alJl2006h . 

3.3. SNR Impacting External Material 

A possibility that has been discussed for SNR J1713 (e.g., 
iBerezhko & Volkll2006Ll2010l) is that the FS is currently im- 
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FIG. 6. — Time history of X-ray emission for the same Model C as 
shown in Fig. \5\ In all panels, the solid black curve shows thermal X-rays 
with Coulomb electron temperature equilibration and the dashed green curve 
shows thermal X-rays with instant equilibration. The solid blue curve is the 
synchrotron continuum. We limit our thermal X-ray calculation to energies 
between 10~ 4 and 10 — 2 MeV, creating the sharp drop-offs beyond these en- 
ergies. The Suzaku observations are shown in red. 

pacting a pre-SN shell or molecular cloud at a radius of 
~ 10 pc. In this case, time scales considerably less than £snr 
become important as the FS moves into the steep density ramp 
of the shell or cloud edge. Since there is a time delay between 
the shock heating of the plasma and the non-equilibrium pro- 
duction of X-ray emission lines, the question naturally arises: 
is it possible for GeV-TeV emission to be produced before 
strong X-ray emission lines are generated? 

In Fig. [5] we show our model C as a time sequence of the 
SNR density profile and the photon emission as the FS runs 
into a density gradient. The parameters are similar to model 
B except we have moved the inner edge of the dense shell to 
~ 9.6pc so that the SNR radius is about lOpc at 1600 yr to 
be consistent with observations. We have also lowered K ep to 
10~ 3 and increased the wind speed to 30 km s _1 from 20 km 
s _1 and made other changes (see Tables[T|[2]and[3]l in the in- 
put parameters in the attempt to obtain a good fit to the X-ray 
and GeV-TeV observations at 1600 yr with pion-decay, from a 
combination of trapped and escaping CR protons, dominating 
the GeV-TeV emission at £snr = 1600 yr. As is clear from 
Figs. |2] and |U the mass concentrated in the dense shell must 
be ^> 100 Mq to make the pion-decay contribution from the 
escaping CRs comparable to IC and we use M s heii = 10 4 M© 
for Fig. 

In Fig. [5] the black dashed curves in the left-hand plots are 
the initial density profile and the red solid curves are the den- 
sity profiles at the ages shown. The right-hand plots show 
the broadband emission at the indicated ages. The important 



properties of Fig. [5] are: 

(i) in order to not overproduce the thermal X-rays, the den- 
sity upstream of the FS at ~ 1600 yr must be < 1 cm~ 3 , forc- 
ing the FS to be ascending the density ramp and not in the 
high-density plateau; 

(ii) the pion-decay emission from the escaping CRs varies 
weakly with isNR since it depends largely on the constant 
amount of external mass in the dense shell; 

(iii) the pion-decay emission from the trapped CRs in- 
creases rapidly as the FS ascends the dense shell because it 
scales ~ rip; 

(iv) as the FS moves into the "molecular cloud," the X- 
ray synchrotron and TeV IC increase less rapidily than the 
pion-decay from the trapped CRs because radiation losses in- 
crease due to the increasing magnetic field in the cloud. This 
is particularly noticeable in the synchrotron between 1400 and 
1600yr; 

(v) between 1600 and 1800yr, the relative intensity of ther- 
mal X-rays increases by more than a factor of ten relative to 
the nonthermal continuum emission; and 

(vi) the distance between the forward and reverse shocks 
drops substantially as the FS runs up the density gradient. 

In Fig. [6] we show the thermal X-rays along with the syn- 
chrotron emission (solid blue curves) in the X-ray band for 
Model C. The solid black histograms are calculated assum- 
ing that electrons immediately behind the shock are cold and 
then heated via Coulomb collisions with the shock-heated 
protons. The dashed green histograms are calculated assum- 
ing the electrons are heated instantly with the protons. While 
the method of electron heating produces important differences 
in the details of the emission lines, all major conclusions 
for broadband models ar e independent of this heating (see 
lEllison etalJ 120071 120 Id for a full discussion of our model 
for electron heating). The calculation of the emission lines 
includes the non-equilibrium effects as the plasma is shock 
heated and compressed and then expands and cools behind 
the FS. For this set of parameters, until ~ 1600 years, the 
emission lines lie mainly at or below the synchrotron contin- 
uum while at 1800yr, the line flux is more than 10 times the 
synchrotron flux at ~ O.lkeV. The variation in line intensity 
between 0.2 and 1 keV over the 400 yr period from 1400 to 
1800yr is greater than an order of magnitude, while the vari- 
ation in synchrotron intensity in this energy range is less than 
a factor of three. 

Also evident in Fig. [6] is that the high-energy synchrotron 
emission decreases with increasing age. We note that 
iPatnaude et al.l (1201 ll) have observed a decline in the nonther- 
mal emission from Cassiopeia A on time scales < 10 yr with a 
1.5%— 2% yr -1 decline in the 4.2-6.0 keV range. Such rapid 
changes in synchrotron emission could result from a number 
of causes as the FS enters a steep density gradient. These 
including the slowing of the FS, an increase in the ambient 
magnetic field strength, as is the case in Fig. [5] and/or the 
damping of magnetic turbulence in the weakly ionized dense 
material. 

3.4. High-Velocity Pre-SN Wind 

Another possibility is that SNR J 17 13 exploded in a 
fast wind with V W i n d ~ 1000km s _1 , and dM/dt = 
10~ 7 Mfn yr -1 , typical of an early B-type or late O-type star 
(e.g.. Ide Jager et ani!988l) . If we accept the constraints that 
the remnant currently has an age of ~ 1600 yr, a radius of 
<10pc, and an explosion energy of i?sN ~ 10 51 erg, it's 
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FIG. 7. — The red curves show the results at tsNR = 1630 yr for a fast 
pre-SN wind where any high-density shell is at a radius beyond 20 pc. The 
model producing the black curves has exactly the same parameters only now 
the dense shell is placed at i? s i lc ii ~ 9.5 pc. In this case, the FS impacts the 
shell well before 1630 yr and is shown 800 yr after the explosion. The arrow 
on the horizontal axis indicates the FS position at 800 yr for Model D. Since 
photon spectra are not calculated for Models D and E, only values important 
for the density profile are given in Tables[2][3] and|4] 

clear that the only possibility is that the FS is now impact- 
ing dense material which stands ~ 10 pc from the explosion 
center. 

To illustrate this, we show in Fig. [7] examples (Mod- 
els D and E), where V w - m d = 500km s _1 , dM/dt = 
10~ 6 M Q yr" 1 , and Af cj ~ 10 M Q . The position of the FS 
at 1600yr will be well beyond lOpc unless the FS is slowed 
by dense material. We note that we have chosen extreme val- 
ues for Fig. 123 a faster V w - m d, a lower dMj dt, and/or a smaller 
M e j, will result in a FS beyond 20 pc at 1600yr. The slower 
speed we use also encompasses the case where the CSM con- 
tains regions of both high and low speed pre-SN winds. 

The black curves in Fig. Q show the case where the shell 
edge is at i? s hcii ~ 9.5 pc, while the red curves show the re- 
sults for the same input parameters except that the dense shell 
is beyond the radius the FS obtains at 4snr = 1630 yr. The 
same arguments used to exclude the FS penetrating the dense 
shell in Fig. [5] apply for a fast wind. The constraints of age, 
distance, and radius are not consistent with the lack of X-ray 
line emission if the FS is interacting with a dense shell. 

4. DISCUSSION AND CONCLUSIONS 

We have generalized our CR-hydro-NEI model of an evolv- 
ing SNR to include the production of radiation in a non- 
homogeneous CSM such as would exist with a pre-SN wind 
or for a SN exploding near a molecular cloud. The model in- 
cludes the production and escape of high-energy CR protons, 
as well as CR electrons and protons that remain trapped within 
the FS. As the escaping CR protons diffuse in the CSM, they 
produce 7-rays by pion-decay and this production is added 
to the synchrotron, inverse-Compton, and pion-decay emis- 
sion from the trapped CRs. As in our previous work (e.g . , 
Ellison et al.ll2007t iPatnaude et all l2009t lEilison et al J 12010) . 
we simultaneously calculate the thermal X-ray line emission 
with the broadband continuum allowing strong constraints to 
be placed on some of the many parameters required to model 
the SNRJ1713. 

Cosmic-ray protons that escape the SNR and interact with 
external material produce 7-rays without producing thermal 
X-rays. If escaping CRs are not considered, homogenenous 
models clearly favor IC from trapped CR electrons as the 



mechanism responsible for the GeV-TeV emission in SNR 
J1713. Including escaping CRs opens the possibility that the 
constraint on ambient density imposed by the lack of X-ray 
line emission in SNR J 17 13 can be satisfied and still have the 
GeV-TeV emission produced predominately by CR protons in 
the external material. 

As in our previous work, we don't calculate the thermal X- 
ray emission at the remnant reverse shock (RS), nor do we 
explicitly include the acceleration of elements heavier than 
protons at the FS. We do include a 10% by number contri- 
bution of helium in the t arget nuclei. The calculations of 
ICaprioli. Blasi. & Amatol (1201 ll) include heavy nuclei in the 
shock modification and show an increase in the pion-decay 
emission by as much as a factor of 2.5 along with changes in 
the shape of the emitted spectrum. However, the enhancement 
of pion-decay emission from heavy elements is not enough to 
modify our main conclusion that leptons dominate the GeV- 
TeV emission, even without considering the thermal emission 
from the RS. Including the RS thermal emission would pro- 
duce an even greater enhancement of IC over pion-decay at 
GeV-TeV energies. 

We have considered four cases. The first is when the exter- 
nal material is a shell from a slow pre-SN wind, in which case 
the external mass is < 100 Mq. The second is when a larger 
external shell of M s heii = 10 4 M Q lies near but outside of the 
FS at tsNR = 1630 yr. The third case examines what happens 
when the FS is interacting with the dense 10 4 M Q shell. The 
fourth case considers a SNR in a fast pre-SN wind interacting 
with a 100 Mq shell. 

Our major conclusion is that it is not possible to obtain a 
satisfactory broadband fit to SNR J1713 dominated by pion- 
decay if the remnant is isolated. An excellent fit can be ob- 
tained (Fig. 123* with IC dominating the GeV-TeV emission, 
confirming that an isolated SNR J1713, whether from a core- 
collapse or thermonuclear SN, produces most of its GeV-TeV 
emission by IC radiation. It is possible, but unlikely, for pion- 
decay from escaping CRs to dominate the GeV-TeV emission 
if this remnant is near or interacting with > 10 4 M of exter- 
nal material. 

If the SNR is near a massive region ~ 10 4 M Q , then 
pion-decay from escaping CRs can, in principle, outshine the 
inverse-Compton from trapped CR electrons at TeV energies. 
However, there are fundamental uncertainties in the proper- 
ties of escaping CRs having to do with the shape of the escap- 
ing CR distribution and their diffusion in the partially ionized 
CSM. The model we use produces a distribution that is too 
narrow to match the combined Fermi-LAT and HESS obser- 
vations with pion-decay alone and a substantial contribution 
from IC is required for a satisfactory fit see Fig.H The high- 
est energy CRs, whether trapped or escaping, will have highly 
anisotropic distributions and substantial work is needed to 
better understand wave generation in such distributions before 
a firm prediction for the shape can be made. 

Our conclusion that the GeV-TeV emission from an iso- 
lated SNR J 17 13 must be dominated by IC is essentially in- 
dependent of the CR diffusion parameters in the CSM since 
escaping CRs contribute an insignificant fraction of the emis- 
sion for M s heii < 100 M©. However, if the escaping CRs 
are interacting with enough mass (> 10 4 Mq), the diffusion 
properties will be important. In that case, the energy de- 
pendence of -Dcsm, as well as the normalization, will in- 
fluence the shape and intensity of the pion-decay emission. 
As the CRs diffuse through the molecular cloud, Z?csm will 
be determined by wave-particle interactions in a dense, par- 



9 



tially ionized medium (e.g jReville et al.l2007l and references 
therein). A detailed analysis of this propagation is beyond the 
scope of this paper. 

A new aspect of our work is the modeling of the remnant 
as the FS impacts a density gradient as shown in Figs. [5] and 
[6] These figures show that the emission, and particularly 
the relative thermal and non-thermal fluxes, will vary sub- 
stantially on short time scales ~ 100 yr depending on how 
long the shock has been interacting with the cloud. We note 
that an important approximation of our CR-hydro-NEl model 
will modify this conclusion somewhat. As the FS runs into a 
steep density gradient, the mass flux crossing the shock will 
increase and, depending on the change in shock parameters, 
more mass will be injected into the DSA mechanism. In real- 
ity, it will take some time for this additional mass to be accel- 
erated to TeV energies but we assume the acceleration occurs 
essentially instantaneously. The X-ray line emission is calcu- 
lated consistently with time so the relative increase in thermal 
over TeV emission shown in Fig. would actually be greater 
than what is shown. 

Significantly, we have found that the core-collapse scenario 
presented here affords a better fit to the broadband SNR J1713 
observati ons than the homog eneous Type la SN model pre- 
sented in Ellis on et al.l (120101) . The lower pre-SN wind mag- 
netic field present in the core-collapse model, yields a lower 
post-shock field, even with strong MFA, allowing electrons to 
be accelerated to higher energies than in the homogeneous 
ISM. These higher energy electrons produce IC emission 
against the CMB consistent with the highest energy HESS 
points. 

The low magnetic field value we find is not necessarily in- 
consiste nt with the rapid time variations observed for SNR 
J1713. lUchiyama et al.1 (120071) observed ~ lyr variations 
in X-ray synchrotron emission and interpreted this as the ra- 
diation loss time scale for TeV electrons in ~ 1 mG fields. 
H owever, other estimates yiel d lower values (see references 
in lEllison & Vladimirovl2008l) . and there is an alternative ex- 
planation fo r rapid time variations that doesn 't require such 
large fields. iBykov. Uvarov. & Ellison! (120081) showed that a 



steeply falling electron distribution in a turbulent magnetic 
field can pro duce intermittent synchr otron emission consis- 
tent with the lUchivama et alJ (120071) observations. In the 
IBykov et al.1 J2008) model, root-mean-square fields of 10's of 
fiG are adequate to explain this time variation and are consis- 
tent with our findings here. 

It is important to note that while we feel the case is all but 
closed for SNR J 17 1 3, we make no claim that all SNRs will 
show IC dominance at GeV-TeV energies. The relative bright- 
ness of pion-decay versus IC depends largely on environmen- 
tal parameters and pion-decay may well dominate in some 
SNRs, particularly if they show strong X-ray line emission. 
It is also possible that, in a particular remnant, high-energy 
emission in some areas will be dominated by IC and in others 
by pion-decay. 

A final important point is that even though IC from rela- 
tivistic CR electrons dominates pion-decay from relativistic 
ions, all of our models, and all nonlinear models we are aware 
of, show the efficient acceleration of protons. Far more en- 
ergy is put into relativistic protons than into electrons by our 
DSA models. Our work does not in any way call into question 
the widely held believe that SNRs are the primary sources of 
CRs, at least up to the so-called knee at 10 15-16 eV. The evi- 
dence that some individual SNRs produce CR ions with high 
efficiency, while indirect is compelling (see ; for example, 
Reynolds & Ellisonl \l992. iHughes et al.l 120001: IWarren et alJ 
2005HHelder et al.ll2009HMorlino & Capriolill201 11) . 
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TABLE 1 

Input Parameters for SN and SNR. 



Model 11 


*SNR 




dM/dt 


^wind 


'-''wind 


rp b 
- 1 wind 




[yr] 


[M Q ] 


[M Q yr" 1 ] 


[kms" 1 ] 




[K] 


A 


1630 


3 


lxlO" 5 


20 


0.03 


10 4 


B 


1630 


3 


lxlO" 5 


20 


0.03 


10 4 


C 


1800 


3 


6xl0~ 6 


30 


0.03 


10 4 


D 


800 


10 


lxlO" 6 


500 


5xl0~ 3 


10 s 


E 


1630 


10 


lxlO" 6 


500 


5xl0 -3 


10 5 



All models use dsNR = 1 kpc, £?sn = 1 x 10 erg, and n = 7. 

The unshocked temperature has very little influence on the solutions as long as it is < 10 6 K. 



TABLE 2 

Input Parameters for DSA and Line Emission. 



ModeP 


£dsa 
[%] 


■^arap 








/norm 


A 


25 


8.5 


0.01 


0.75 


1 


0.75 


B 


25 


8.5 


0.01 


0.75 


1 


0.65 


C 


50 


10 


10" 3 


0.75 


0or 1 


0.9 


D 


25 












E 


25 













All models use / s k = 0.1, N C ut = 30, and a solar elemental composition. 



TABLE 3 

Input Parameters for CSM and escaping CR diffusion. 



Model' 1 


D CSM,0 


Acsm.o 


"shell 


Afshcll 


Bshell 


-Bshcll 




[cm 2 s ] 


[pc] 


[cm" 3 ] 


[Mq] 


[pc] 


[AtG] 


A 


1 x 10 26 


3.2 xlO" 3 


1 


LOO 


12 


3 


B 


1 x 10 26 


3.2 xlO" 3 


100 


10 4 


12 


3 


C 


lxlO 27 


0.032 


10 


10 4 


9.6 


1 


D 






1 


100 


9.5 




E 






1 


100 


> 20 





All models use a rg = /3 n = 0.5, and n un ; = 0.01 cm 3 . 



TABLE 4 

Output Values at end of simulation. 



Model 


rtot 


n p a 




B 2 b 






-/Vis wept 




6esc 


Xinj 






[cm -3 ] 


L«G] 


I/^G] 


[pc] 


[kms" 1 ] 


[M ] 


[E SN ] 


[E S n] 




A 


4.6 


0.013 


0.22 


10 


S.9 


4200 


4.2 


0.15 


0.013 


3.87 


B 


4.6 


0.06 


0.23 


10 


S.9 


4200 


4.2 


0.15 


0.013 


3.87 


C 


5.6 


2.0 


0.7 


75 


9.8 


850 


4.3 


0.17 


0.03 


4.0 


D 


4.6 


0.01 






9.7 


2500 


0.025 


4xl0" 3 


3xl0- 4 


3.9 


E 


4.6 


1.5xl0~ 5 






19 


9000 


0.04 


6xl0 -3 


5xl0" 4 


3.7 



'Value just upstream of the FS at the end of the simulation. 

'Magnetic field just downstream from the FS at the end of the simulation. 



